X  AD-A092  475  SCIENCE  APPLICATIONS  INC  LA  JOLLA  CA  F/6  13/10 

TIHE-OONAIN  METHOD  FOR  COMPUTING  FORCES  AND  MOMENTS  ACTING  ON  T— ETC(U) 
SEP  80  R  B  CHAPMAN  NOOO14-70-C-O111 

UNCLASSIFIED  SAI-462-80-560-LJ  NL 


SAI  REPORT  NUMBER 
SAI-462-80-560-LJ  J 


£  -  TIME-DOMAIN  METHOD  FOR  COMPUTING 
FORCES  AND  MOMENTS  ACTING  ON 
THREE  DIMENSIONAL  SURFACE-PIERCING 
SHIP  HULLS  WITH  FORWARD  SPEED  , 


by 

Richard  B. 


apman 


Science  Applications,  Inc. 

1200  Prospect  Street 
La  Jolla,  California  92038 


/ 


Sepl 


Final  ^ep«Pt .  15  Dec* 


177—  3<}>  Sepl 


APPROVED  FOR  PUBLIC  RELEASE 
DISTRIBUTION  UNLIMITED 


Prepared  for 

DAVID  W.  TAYLOR  NAVAL  SHIP  R&D  CENTER 
Bethesda,  MD  20084 

OFFICE  OF  NAVAL  RESEARCH 
800  N.  Quincy  Street 
Arlington,  VA  22217 


.  -  v 

Jf  \  " 


80. 


UNCLASSIFIED 


security  classification  or  this  pa ge  /whah  o««  fn»*r»« 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

'■  REPORT  NUMICA  j.  OOVT  ACCESSION  NO. 

SAI-462-80-560-LJ  ^  0' ft 6^  V/^ 

»•  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  tm*  Minim) 

TIME-DOMAIN  METHOD  FOR  COMPUTING  FORCES 
AND  MOMENTS  ACTING  ON  THREE  DIMENSIONAL 
SURFACE-PIERCING  SHIP  HULLS  WITH  FORWARD 
speed 

i.  tyre  of  report  a  perioo  covered 

Final:  15  Dec  1977  - 

30  Sep  1980 

4.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHOR/ *| 

R.  B.  Chapman 

•-  contract  or  grant  NUMBER/*; 

N00014-78-C-0111 

*■  performing  organization  name  ano  address 

Science  Applications,  Inc. 

1200  Prospect  Street 

La  Jolla,  CA  92038 

10.  PROGRAM  Clement,  project,  task 
AREA  4  WORK  UNIT  NUMBERS 

122301 

1 1.  CONTROLLING  office  name  ano  aooress 

David  W.  Taylor  Naval  Ship  R&D  Center 
Bethesda,  MD  20084  (Code  1505) 

It.  report  DATE 

September  1980 

IS.  NUMBER  OF  PAGES 

34 

14  monitoring  AGENCY  NAME  4  AOORCSS/II  millmrmmi  (rmm  CmnrmlUnt  Olllem) 

Office  of  Naval  Research 

800  Quincy  Street 

Arlington,  VA  22217 

IS.  SECURITY  CLASS,  (ml  Hum  impmtl) 

Unclassified 

lia.  declassification  OOWNGRAOING 
schedule 

<4.  DISTRIBUTION  STATEMENT  (ml  INI*  StMrII 

APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 

IT  DISTRIBUTION  statement  (ml  IA*  aAairac I  miwM  Mi  Urn* 4  TO.  Il  AIIhmi  Irmm  Hmmmrf) 

■  •  surrlementary  notes 

Sponsored  by  the  Naval  Sea  Systems  Command,  General  Hydromechanics 
Research  Program  administered  by  David  W.  Taylor  Naval  Ship  Re¬ 
search  and  Development  Center,  Code  1505,  Bethesda,  MD  20084 

••  NET  ROROS  (Cmtnmum  mm  rmmmrtm  mtmm  II  moimt  in  immmilly  4r  414*4  mmmrnmr) 

Ship  motions,  Seakeeping,  Hydrodynamics,  Numerical  Fluid  Dynamics 

1 

\ 

M  abstract  rCMiMM  «  rmwmmmm  NR  II  MiHMir  m4  Hmmlltr  *r  W«»  ■ m*n| 

‘a  time-domain  simulation  method  for  computing  forces 
and  moments  acting  on  an  arbitrary  surface-piercing  three-dimen¬ 
sional  ship  hull  is  presented.  Arbitrary  motions  can  be  pre¬ 
scribed  and  are  assumed  to  be  sufficiently  small  so  that  the 
linearized  method  is  valid.  Forward  speed  effects  are  included 
under  the  assumption  that  the  disturbance  generated  by  forward 
motion  is  also  small  and  interactions  with  the  flow  generated  by 

DO  .'2Tn  1471  edition  of  •  rov  ••  it  oiwkiTt  UNCLASSIFIED 


mcdoitt  classification  >»  tw>  »m  iK> WSSVmmm 


_ UNCLASSIFIED _ 

*CCUN|TV  CLASSIFICATION  OF  THIS  FAOEfWMw  Df  Kitlf*) 

"  A  | 

the  prescribed  motions  are  of  second  order.  The  hull  is  repre¬ 
sented  by  a  set  of  quadrilateral  surface  panels  in  a  body-fixed 
system,  while  the  free  surface  is  represented  by  its  spectral 
coordinates  in  a  space-fixed  rectangular  system.  Small  time  steps 
are  used  to  advance  the  flow.  This  time-domain  computation  can 
also  be  applied  to  compute  linearized  seakeeping  responses  in 
the  frequency  domain.  For  each  degree  of  freedom,  a  unit  step 
function  is  applied  to  the  velocity.  The  resulting  impulsive 
forcds  and  moments  and  the  computed  time  histories  of  these  forces 
and  moments  after  the  step  in  velocity  can  then  be  combined  to 
compute  the  added  mass  and  damping  forces  and  moments  in  the 
frequency  domain.  Wave-excited  forces  can  be  computed  by  a  simi¬ 
lar  method.  As  a  check  on  the  computer  program,  the  added  mass 
and  damping  of  a  semi-submerged  sphere  oscillating  in  heave  were 
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lytic  mass  and  damping  coefficients  for  a  fairly  coarse  represen¬ 
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ABSTRACT 


A  time-domain  simulation  method  for  computing  forces 
and  moments  acting  on  an  arbitrary  surface-piercing  three-dimen¬ 
sional  ship  hull  is  presented.  Arbitrary  motions  can  be  prescribed 
and  are  assumed  to  be  sufficiently  small  so  that  the  linearized 
method  is  valid.  Forward  speed  effects  are  included  under  the 
assumption  that  the  disturbance  generated  by  forward  motion  is 
also  small  and  interactions  with  the  flow  generated  by  the  pre¬ 
scribed  motions  are  of  second  order.  The  hull  is  represented  by 
a  set  of  quadrilateral  surface  panels  in  a  body-fixed  system,  while 
the  free  surface  is  represented  by  its  spectral  coordinates  in  a 
space-fixed  rectangular  system.  Small  time  steps  are  used  to 
advance  the  flow.  This  time-domain  computation  can  also  be 
applied  to  compute  linearized  seakeeping  responses  in  the  frequency 
domain.  For  each  degree  of  freedom,  a  unit  step  function  is 
applied  to  the  velocity.  The  resulting  impulsive  forces  and 
moments  and  the  computed  time  histories  of  these  forces  and  moments 
after  the  step  in  velocity  can  then  be  combined  to  compute  the 
added  mass  and  damping  forces  and  moments  in  the  frequency  domain. 
Wave-excited  forces  can  be  computed  by  a  similar  method.  As  a 
check  on  the  computer  program,  the  added  mass  and  damping  of  a 
semi-submerged  sphere  oscillating  in  heave  were  computed  numeri¬ 
cally.  The  results  show  good  agreement  with  analytic  mass  and 
damping  coefficients  for  a  fairly  coarse  representation  of  the 
hull  and  free  surface. 


INTRODUCTION 


This  report  describes  a  computational  method  designed  to 
simulate  arbitrary  linearized  motions  of  a  surface-piercing  three- 
dimensional  ship  hull  in  the  time  domain.  The  hull  is  specified 
by  a  set  of  quadralateral  surface  panels  describing  the  hull  por¬ 
tion  directly  beneath  the  static  waterline.  The  numerical  method 
calculates  the  forces  and  moments  generated  by  arbitrary  motions 
of  the  body  and/or  the  action  of  an  ambient  wave  field  acting  on 
the  body  under  the  assumption  that  the  flow  velocities  can  be 
linearized. 

The  method  may  be  classified  as  a  time-domain  computation 
analogous  to  the  frequency-domain  methods  developed  by  Chang  (1) 
for  three-dimensional  ship  hulls.  The  time-domain  method  is  the 
more  fundamentally  numerical  approach  to  this  class  of  problems. 

The  time-domain  method  does  not  use  Green's  functions  representing 
solutions  of  the  free-surface  equations  for  singularities  follow¬ 
ing  prescribed  conditions  (i.e.,  moving  with  uniform  speed  and 
oscillating  as  in  reference  [1]).  Thus  the  time  domain  method 
provides  the  flexibility  of  arbitrary  motions.  It  can  easily  be 
generalized  to  include  a  time-dependent  hull  shape  or  other 
generalized  features. 

Limitations  on  computer  times  and  storage  do  not,  however, 
allow  this  "brute  force"  type  method  to  serve  as  a  practical  tool 
in  three  dimensions  at  the  present  time.  With  continued  improve¬ 
ments  in  computer  technology,  however,  methods  of  this  nature  may, 
in  the  near  future,  be  explored  for  practical  applications.  It  is 
hoped  that  the  techniques  described  here  will  provide  a  framework 
for  later  work  of  this  nature. 

A  basic  objective  in  the  development  of  the  computational 
method  was  to  establish  the  method  in  as  simple  a  form  as  possible. 
The  resulting  numerical  technique  described  in  following  sections 
of  this  report,  is  a  generalization  of  the  two-dimensional  method 
applied  in  reference  [ 2] ,  expanded  to  three  dimensions  with  some 


improvements  in  the  numerical  method.  Another  objective  was  to 
examine  the  possible  problem  areas  such  as  the  free  surface  area 
interior  to  the  hull  and  pressure  computations.  These  areas  are 
also  discussed  in  this  report.  Direct  application  to  linearized 
seakeeping  computations  is  discussed.  Finally  a  sample  problem  - 
heave  oscillations  of  a  semisubmerged  sphere  -  is  presented. 

Problem  Formulation 

Consider  a  coordinate  system  (x,  y,  z,  t)  with  z  positive 
downward  and  z  =  0  at  the  static  waterline.  Far  from  the  origin 
the  flow  is  at  rest.  This  space-fixed  coordinate  system  is  the 
basis  for  describing  the  fluid  flow  in  general  and  the  free  surface 
in  particular.  The  velocity  potential  <j)(x,  t)  satisfies  Laplaces 
equation  throughout  the  fluid. 


<j>  =  0 


z  >  0 


The  free  surface  elevation  field  n(x,  y,  t)  and  its  time  derivative 
n(x,  y,  t)  are  defined  over  the  entire  plane  z  *  0.  Their  values 
at  any  time  provide  sufficient  initial  conditions  to  define  the 
free  surface  problem.  Exterior  to  the  intersection  of  the  hull 
surface  the  elevation  field  satisfies  the  dynamic  and  kinematic 
free  surface  conditions 


3<fr  _  3n 

3z  ~  3t 

On  the  free  surfaces  interior  to  the  hull  any  condition  could  be 
applied.  In  this  case,  however,  we  assume  that  the  above  linearized 
equations  apply  over  the  entire  z  =  0  plane.  An  arbitrary  surface¬ 
piercing  hull  is  defined  by  the  surface, 


F(x,  y,  z,  t) 


z  >  0  , 


l 


2 


r 

with  the  exciting  normal  velocity  vn  induced  by  body  motions  or 
interactions  with  the  ambient  wave  field  prescribed  on  the  surface 
of  the  hull, 

vQ(x,  t)  =  v*(x,  t)  *n,  x  £  shull  • 

(In  practice  the  normal  acceleration  is  specified  rather  than  the 
velocity. )  The  disturbance  generated  by  the  body  motion  is  divided 
into  two  parts 


*  "  ^BD  +  4>FS 


where  t})^  represents  the  instantaneous  effect  of  the  body,  while 
<()pg  represents  the  free  surface  disturbance  generated  by  the  body 
over  all  previous  motions  or  interactions  with  the  ambient  wave 
field.  In  particular  the  boundary  conditions  at  the  free  surface 
(z=0)  are 


and  on  the  body. 


^BD 

34>pg 

at 


o  z  =  o 
-gn  z  =  o 


—  (<f) 
an 


BD 


<f>FS)  =  v**n  x  e  S. 
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It  is  assumed  that  body  components  velocities  are  small,  i.e., 


vbd(x)  =  0  (e) 

and 

vpg(x)  =  0  (e)  • 

2 

Any  terms  of  order  £  are  dropped.  For  example,  this  analysis 
applies  to  the  problem  of  ship  motion  in  a  seaway  at  finite  speed 
if  it  can  be  assumed  that  the  flow  field  induced  by  forward  motion 
alone  and  the  flow  field  induced  by  the  seaway  acting  on  the  ship 
are  both  small  and  do  not  interact.  In  this  case  of  finite 
forward  speed,  the  velocity  at  any  point  on  the  hull  x  is 
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vs(x)  =  U  ex  +  0  (e) 

where  U  is  the  forward  speed.  The  disturbance  and  resulting  forces 
generated  by  forward  speed  alone  are  steady  (for  steady  u)  and  are 
not  addressed  here  (although  this  method  could  be  easily  applied). 
It  is  assumed  that  the  velocity  field  generated  by  forward  speed 
alone  is  of  order  e. 


Outline  of  Numerical  Method 

Starting  from  a  condition  of  rest,  the  simulation  is 
achieved  by  a  series  of  small  time  increments.  At  each  time  step, 
three  fundamental  sets  of  variables  -  which  essentially  define  the 
flow  at  the  instant  -  are  computed; 

(a)  The  time  rate  of  change  of  the  strengths  of  the  hull 

source  panels  in  a  hull-fixed  system  are  computed  from 
the  hull  acceleration  boundary  condition.  This  deter¬ 
mines  the  time-derivative  of  the  body-induced  component 
of  the  flow,  (as  well  as  <|>BD  through  integration). 

(b)  The  free  surface  elevation  field  n(x,  y)  and  its  time 

derivative  n(x,  y)  are  computed  by  advancing  the  free 
surface  by  a  small  increment  in  time  and  adding  the  wave 
elevations  generated  by  the  hull  source  panels  over  that 
increment.  This  defines  the  free  surface  component  of 
the  flow,  and  its  time  derivative  in  the  space-fixed 

system. 

(c)  Pressures  acting  on  the  hull  are  computed  separately  for 
the  free  surface  and  body  components:  pBg,  PBD-  Both 
are  computed  at  hull  panel  center  points  and  added  to 
yield  total  forces  and  moments. 
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Computation  of  Panel  Source  Strengths 


The  hull  is  represented  numerically  by  a  set  of  Ng 
quadralateral  panels.  Each  panel  is  specified  by  the  coordinates 
of  the  four  corner  points,  a  panel  center  point,  and  a  normal 
vector.  The  fluid  velocity  at  each  center  point  in  the  normal 
direction  may  be  written  as 

^BD.*Si  +  ^FSi*“i  =  V(t)  *“i  i  =  1,  2,  ...  Ng 

where  v^*(t)  is  a  prescribed  forcing  function  depending  on  the 
motion  of  the  hull  and  the  ambient  wave  field.  Since  the  free 
surface  induced  velocity  v^g  is  assumed  to  be  available  from  the 
free  surface  representation,  this  equation  fixes  v™  •  n.  at  the 

OJJ  1 

center  point  of  each  hull  panel. 


A  simple  source  and  its  image  is  assumed  to  be  distri¬ 
buted  over  each  hull  panel 

G(x,  y,  z,  x",  y " ,  z")  =  ((x-x')2  +  (y-y')2  +  (z-z' )2)-1/2 

-  ( ( x— x ^ ) 2  +  ( y-y ')2  +  (z+z")2)_1//2 


with 


4>„n(x,  y,  z) 


■  A 

1=1  o 


z,  x",  y' ,  z" )dS"  , 


where  is  the  strength  on  the  i  source  panel  and  S7  represents 
the  panel  surface. 

Consider  a  simple  source  of  unit  strength  acting  over  a 
plane  surfaces', 


»  =  #ds' 


X  -  X 


Then  it  may  be  shown  analytically  that  the  velocity  normal  to  the 
plane  is  identical  to  the  solid  angle  occupied  by  the  panel  surface 
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and  the  velocity  tangential  to  the  plane  surface  can  be  written  as 
a  line  integral  around  the  boundary  of  the  panel.  Numerical  appli¬ 
cation  of  these  results  produces  an  array  E. .  giving  the  normal 
velocity  component  at  the  center  point  of  panel  i  induced  by  a 
uniform  source  distributed  over  panel  j  and  its  image, 

NT 

Vr 


B 


'BD.  *  ni  “  £  E 
i  d=l 


.  .  a  . 
iJ  J 


i  =  1,  ...  nb. 


This  provides  a  set  of  linearized  equations  which  can  be  used  to 

determine  the  source  strengths  a..  In  practice,  it  was  found  that 

J 

better  estimates  for  body-induced  pressures  could  be  obtained  by 
using  the  local  acceleration  rather  than  the  velocity  directly. 

The  time  derivative  relative  to  a  hull  fixed  coordinate  system 
is  computed  at  the  beginning  of  each  time  step  from 

Nr 


J*!  Eij  °j  aBD 


ni  ’  ai*(t>  -  aPSl  '  °i  -  3  &  vFSi  •  \ 


where  aFg  and  vFg  are  the  free-surface  induced  acceleration  and 
velocity  ^in  space-fixed  coordinates)  at  point  i  and  ai*(t)  is 
the  prescribed  normal  acceleration  in  hull-fixed  coordinates. 

For  example,  if  the  body  has  surge  acceleration  a^Ct),  then 
ai*(t)  =  a^Ct)  •  nxl*  The  source  strengths  at  the  center  can  be 
computed  directly  from  the  local  velocities  if  desired.  In  prac¬ 
tice  it  was  found  that  they  could  be  computed  adequately  by  inte¬ 
grating  their  derivatives  in  time. 


Representation  of  Wave  Elevation  Field 

The  representation  of  the  wave  elevation  field  and  its 
time  derivative  by  their  spectral  components  allows  the  linearized 
free  surface  boundary  conditions  to  be  applied  easily  when  adapting 
the  free  surface  component  of  the  flow.  It  also  aids  in  computing 
free  surface  induced  pressures  and  velocities  on  the  hull.  An 
arbitrary  real  continuous  function  n(x,  y)  defined  over  the  z  *  0 
plane  can  be  written  in  the  form 

00  00 

n(x,  y)  =  J dkx  J ' dky  A(kx,  ky)  ei^kx  *  x  +  kY  '  30 
0  —00 

where  A(kx,  ky)  is  complex  function  and  the  elevation  n(x,  y)  is 
taken  to  be  the  real  part  of  the  expression  on  the  right  hand  side-. 
Numerically  n(x,  y)  can  be  represented  by 


n(x,  y) 


NKX  NKY  T 

’  4kxn  m?,  4kym  [ 

n-1  m=l  u 


A  (t)  e 
nm 


i(kx  •  x  +  ky  •  y) 


+  A*  (t)  e 
nm 


i(kx  •  x  -  kyr 


•  y)j 


Similarly  if  n(x,  y)  is  defined  as  the  time  derivative  of  n(x,  y) 
due  to  free  surface  effects  alone,  that  is  the  rate  of  change  of 
n(x,  y)  in  the  absence  of  the  hull,  it  may  be  represented  numeri¬ 
cally  by 


n(x,  y)  = 


NKX  NKY  f 

*  4kxn  4kym  L 

n=l  m=l 


B  ( t )  a  e 
nm  nm 


+  B*  ( t )  a  e 
nm  nm 


i(kxn  •  x  +  kym  •  y) 


i(kxn  •  x  -  kym 


•  y)j 


where  a  >  0  is  defined  by 
nm 


2  ,,  2  .  .  2.1/2 

=  S(kx„  +  ky) 


Thus  in  the  absence  of  the  hull  the  time  derivative  of  A  is 

nm 

given  by 

4-  Aft)  =  B  ( t )  cnm 

dt  nm  nm  nm 


-jr  A*  (t)  =  B*  (t)  o 
dt  nm  '  nm  nm 


From  the  kinematic  free  surface  condition  for  the  free  surface 
component  alone , 


3 r|  _  ^FS 


,2  =  0 


a  special  representation  for  the  free  surface  potential  may  be 
derived 


where 


2)  3  2  Akx 

£ 

Aky  B  ( t  )a  k-3- 

Jm  1  nmv  '  nm  nm 

n=l 

m=l 

w 

+  B*  (t) 

a 

i(kx  *  x  -  ky 
k  1  e  n  m 

nm' 

nm 

nm 

knm 

,2,2 
=  kx„  +  ky 
n  •'m 

i(kxQ  •  x  +  kym  •  y) 


•  y)J 


-k  •  z 
nm 


Now  the  dynamic  free  surface  condition, 

tF  -  -gn  *  -  0, 

can  be  applied  to  yield  the  remaining  equation  to  be  used  when 
updating  the  free  surface  in  time, 

4  B  ( t )  =  -A  ( t )  a 
at  nm  nm  nm 


*nm 
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The  free  surface  induced  pressure  and  normal  component  may  be 

computed  from  A _  and  A*  directly 

r  nm  nm  J 


NKX  NKY  (ikx  •  x  -  k  •  z) 

=  pg  £  Akyn  L  Akym  ®  n  nm 

n=l  n  m=l  m 


,(t)  e 


ikym  •  y 


-iky 

+  A*  (t)  e  m 
nm 


NKX  NKY  (ikx  •  x  -  k  •  z) 

n  -  -g  £  Akx  £  AkYm  e  n 
n=l  n  m=l  m 


,(t)  •  e 


iky^  •  y 


(ikyn  •  nx  +  ikym  •  ny  -  ^  •  nz> 


-iky  •  y 

+  A*  (t)  •  e  •  (ikx  •  n  -  iky  •  y  -  k  • 

nm  n  x  m  nm 


v] 


The  distribution  of  values  of  kx  and  ky  where  A  .  A* 

n  '  m  nm  nm 

and  are  specified  is  critical.  For  practical  computing 

nm  nm 

times,  the  number  of  points  should  be  kept  to  an  absolute  minimum. 

On  the  other  hand,  the  free  surface  representation  must  be  valid 

(1)  Over  the  hull  surface  throughout  the  flow  in  the  immediate 
neighborhood  of  the  hull  surface  and, 

(2)  Over  *-ue  maximum  anticipated  time  interval  of  the 
simulation . 


The  hull  and  the  free  surface  flow  is  assumed  to  be 
located  in  the  general  neighborhood  of  the  origin  of  the  (x,  y) 
plane.  Suppose  the  free-surface  component  of  the  flow  were  periodic 
in  space  with  periodic  lengths  2it  •  Lx  and  2tt  •  Ly  in  the  x  and  y 
directions  respectively.  Then  the  continuum  of  kx  and  ky  are 
replaced  by  the  discrete  values  - 


Similarly,  if  the  free  surface  component  of  flow  were  periodic 
in  time  with  a  period  of  2ttT,  then 


.  2  ^  ,  2,1/2  2.  _2 

k  =  ( kx  +  ky  )  '  =  n  /gT 


These  results  are  applied  in  obtaining  a  kx,  ky  distribution 
satisfying  the  following  requirements.  The  free  surface  flow 
representation  must  be  valid. 

(1)  Over  a  region  extending  L  in  the  x  direction  and  L  in 

x  y 

the  y  direction  from  the  origin  of  the  (x,  y)  plane, 

(2)  Over  features  with  characteristic  lengths  of  &x  and  l 
or  larger 

(3)  Over  time  intervals  of  length  T  or  shorter 


without  significant  influence  from  the  discretization  of  kx  and  ky. 

Condition  (1)  imposes  maximum  step  sizes  for  kxQ  and  kym 


kx  . , 

-  kx  < 

L"1 

n+1 

n  - 

X 

-1 

^m+l 

-  kymi 

V 

Condition  (2)  sets  the  minimum  upper  bounds  for  kxQ  and 


kx 

ky 


> 

A'1 

NKX  - 

X 

> 

A'1 

MKY  - 

y 

Condition  (3)  cannot  be  applied  directly  for  independent 
values  of  kxn  and  kym<  It  can  be  adequately  replaced  by  a  pair  of 
similar  conditions, 

/kxn+i  -  /kxn  £  (T  •  /g)-1 

and 

^  £  <T’ 


These  final  conditions  usually  affect  only  the  few  lowest  values  of 
kx  and  ky,  depending  on  the  values  of  the  other  parameters. 
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Advancement  of  the  Free  Surface  in  Time 

With  each  time  step,  the  free  surface  portion  of  the  flow 

as  represented  by  A^d),  A*m(  t )  and  B*m(t)  h(x,  y)) 

must  be  advanced.  The  new  values  are  those  which  would  result  if 

the  hull  were  not  present  over  the  time  interval  (from  t  to  t  +  At) 

n  n 

added  to  estimates  for  the  free  surface  flow  generated  by  the  body 
(i.e.,  by  the  source  panels  and  their  images)  over  the  small  time 
interval.  The  algorthms  used  to  update  the  coefficients  are 


A  (t  +  At) 
nm  n 


BODY 

A  ( t  )  cosa  At  +  B  ( t  )  sina  At  +  AA 
nm  n  nm  nm  n  nm  nm 


and 


Bit  +  At) 
nm  n 


RODY 

B  (t  )  cosc  At  -  A  (t  )  sina  At  +  AB° 
nm  n  nm  nm  n  nm  nm 


(with  identical  algorthms  for  updating  A*m  and  B^).  Since  the 
body-induced  potential  is  defined  as  zero  on  the  free  surface, 

<J>BD(x,  y,  z)  =  0,  z  =  0, 

the  influence  of  the  body  source  panels  and  their  images  on  Bnm, 

B*  over  this  short  time  period  are  second  order  in  At 
nm 

iBB0DY,  abbody*  ;3ut>2 

nm  ’  nm  ^ 


The  free  surface  elevation  induced  by  the  body  source 
panels  and  their  images  over  a  small  time  interval  At  may  be  com¬ 
puted  by  summing  the  contributions  of  each  panel  and  its  image. 

A  panel  with  uniform  density  a^^  acting  over  a  surface  of  area  A^^ 
is  replaced  by  a  single  point  source  with  strength 

si<t>  *  WV  *  <t-tn)  5(V> 

located  at  the  panel  center  point  (x^  yi(  z i)  as  defined  at  the 
midpoint  of  the  interval,  t  =  t  +  1/2  At. 
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The  vertical  velocity  at  z  *  0  induced  by  these  source 
points  and  their  images  may  be  written  as 


.  £ 

jLl 


si(t) 


i=l  (x-xi)  +  ( y-y± )  + 


2  +  2,3/2  • 


This  can  be  expressed  in  integral  form  as 

2  /*  n  B  r  ikx(x-x.)  r 

An0(x,y)  *  ~  I  dt  £  si(t)/  dkx  e  /  dky  cos  ky(y-yi) 

tn  0 

2  2 

-(kx+ky )z . 

•  e  1 

BODY  BODY* 

Thus,  to  second  order  in  At.AA  _  and  AA  can  be  written  as 

nm  nm 


nr\Tw  **  B  ,  -i(kx*x.  +ky*y.)  -k  z, 

A.B0DY=lOAt  Vo/'*  .1  A*\  «  n  1  ra  l  o  nm  i 

AA  “2  — -  2->  ( t:  +  *  At)  e  e 

nm  n  i»i  i  n  2 


aDd  £  *  1  4t)  e-i(,Iyn-Xi  *  “Vyl>  .**»-  »i 

nm  tt  in  2 

Also  in  second  order  in  At,  the  changes  in  B  induced  by  the  body 

nm 


“>rT  -  -  ^(t) 


onrtv  /•a*\2  B  -i(kx*x.+ky*y.)  -k  z.  a 

ABB0DY  «  -  V  s,(t  )  e  n  1  ym  yi  e  nm  i  nm 

nm  tt  i  n 


with  a  similar  expression  for  AB 


BODY* 
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Computation  of  the  Body-Induced  Pressures 

The  total  force  acting  on  the  hull  surface  is  determined 
from  the  pressures  computed  at  the  center  points  of  the  hull  panels. 
These  pressures  are  assumed  to  act  uniformly  over  the  panel  sur¬ 
face.  From  Bernoulli's  equation,  the  pressure  at  a  point  fixed  in 
space  is 

p  »  _p  |1  +  °  ( e2 ) 

where  the  time  derivative  is  in  a  space-fixed  coordinate  system 
and  the  velocity  field  relative  to  this  space-fixed  system  is 
assumed  to  be  of  order  e.  The  pressure  may  be  written  as  the  sum 
of  free-surface  and  body-induced  components 


where 


P  =  PFS  +  PBD 


PFS  =  -p 


3$FS 

at 


and 

d<^BD  -*■  * 

PBD  =  “d T  +  P  Vs*  7  ^BD 

where  vg  is  the  velocity  of  a  point  on  the  hull  surface.  If  the 
body  moves  with  translational  velocity  Vg  and  rotational  velocity 
Ug  then 


v  = 

s 


B 


Wg  X  X 


The  term 


d<(> 


BD 


dt 


represents  the 


potential  at  a  point  fixed  on  the  hull 
In  most  cases  of  interest 


time  derivative  of  the 
(i.e.  at  a  panel  center). 


vg  *  u  ex  +  0  (e) 
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where  U  is  the  forward  speed  and  e  is  the  unit  vector  in  the  x 
direction.  Unless  otherwise  specified  the  body-induced  pressure 
is  assumed  to  be 


=  -p 


+  pU 


The  time  derivative  — at  the  center  point  of  panel  i  can  be 
calculated  from  the  previously  computed  time  derivative  of  the 
source  strengths  in  a  hull-fixed  system, 


“♦bd/V 


*  E  p.,  “.(V 

jTi  j  n 


where  P.  .  represents  the  potential  at  the  center  of  panel  i  induced 

by  a  uniform  source  density  of  unit  strength  acting  over  panel  j. 

Unfortunately,  elements  of  the  matrix  P.  .  giving  the  potential  are 

i  J 

not  as  easy  to  compute  numerically  as  the  matrix  E.  .  giving  the 

J 

velocity.  At  present  integrals  of  the  form 


i  \  |x.-x|  |x.-x  I  / 


where  x  is  the  image  of  point  x  on  surface  S.,  are  estimated  by 

J 

dividing  the  surface  S .  into  many  TSmall  elements  and  integrating 

J 

numerically  with  the  integrand  evaluated  at  the  center  of  each 
subelement.  This  'brute  force'  method  can  consume  significant 
computing  time  if  resonable  accuracy  is  desired.  These  elements 
need  be  computed  only  once,  however,  if  the  hull  does  not  change 
over  the  time  period  of  simulation. 

The  forces  and  moments  acting  on  the  body  in  a  hull- 
fixed  coordinate  system  may  be  written  in  the  form, 


-  £  "mi  Ai  Pi 


m  =  1 
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where  Fm  is  a  generalized  force  corresponding  to  the  m 


th 


degree  of 


freedom,  n  .  is  the  generalized  normal,  A.  is  the  panel  area  and  p. 

in -L  . .  l  l 

is  the  pressure  on  the  i  panel.  From  this  form,  the  force  F' 

d4>BD  m 

resulting  from  the  — g-  -  term  can  be  written  as 

Ns 

Fm  “  p  ^  nmi  Ai  dt  <pBDi 


NB  nr 

S  Ai  £  pu 


nb 

=  p  £  pp  •  o .  . 
jtl  — 1 J 


The  array  of  (6  x  Ng)  elements  can  be  evaluated  once  at  the 

beginning  of  the  computation,  (t  =  t  ). 

n 

Similarly  the  methods  used  to  calculate  the  matrix 

can  also  give  elements  of  a  matrix  X. .  such  that 

^  J 


N. 


Vx  ^BD . 


B 


(t  ) 


E  ctn)  • 

jti  ^  J  n 


The  forces  resulting  from  the  second  term  of  the  eruation  for  p 


BD 


are 


N. 


N. 


Fm 


'°Uj?i  -  Ai  j?iXiJ  ^  (V 

nb 

-pU  E  PX.  a. 
jtl  — J 


where  PX  .  is  a  (6  x  NB)  matrix  similar  to  PP  .  . 
— raj  — mj 
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Alternative  Techniques 

Certain  subelements  of  the  computation  can  be  represented 
by  alternative  forms.  The  computational  method  described  above 
was  chosen  from  the  following  set  of  possibilities 

(1)  Either  source  or  dipole  panels  can  be  used  to  represent 
the  hull  surface.  If  dipoles  are  used,  then  the  normal 
velocity  on  the  image  panel  is  the  negative  of  the  vel¬ 
ocity  on  the  panel  itself  so  that  the  net  mass  inflow 
into  the  region  enclosed  by  the  hull  and  its  image  is 
zero.  The  major  numerical  advantage  of  dipole  panels 
is  that  the  potential  induced  by  a  uniform  dipole 
strength  of  unit  magnitude  is  easily  computed  at  any 
arbitrary  point  in  space  (it  is  simply  the  solid  angle 
enclosed  by  the  panel  relative  to  that  point).  This 
allows  one  to  compute  the  pressures  induced  by  the  time 
derivative  of  the  body-induced  potential  much  more  easily 
numerically  than  if  source  panels  were  used.  The  primary 
disadvantage  of  dipole  panels  is  that  the  free  surface 
disturbance  generated  by  them  is  discontinuous  along 
points  where  the  hull  intersects  the  free  surface.  That 
is,  there  is  a  jump  in  elevation  between  the  interior 
and  exterior  regions.  This  requires  much  finer  resola- 
tion  in  the  representation  of  the  free-surface  and,  there¬ 
fore,  more  points  in  k-space  for  the  spectral  distribu¬ 
tion  of  the  free-surface. 


(2)  Alternative  conditions  can  be  applied  to  the  free  surface 
interior  to  the  hull.  For  example  instead  of  using  the 
linearized  free-surface  condition  after  the  entire  free 
surface,  the  surface  elevation  interior  to  the  hull  can 
be  constrained  to  be  uniformly  zero  by  placing  surface 
panels  over  this  region.  One  possible  advantage  of 
"adding  a  lid"  is  that  it  suppresses  standing  waves  which 
otherwide  are  excited  in  the  hull  interior.  The  major 
disadvantage  of  a  "lid"  is  the  extra  computing  time 


(3) 


required.  In  theory  these  internal  standing  waves  pro¬ 
duce  no  net  pressure  on  the  (exterior)  hull  surface. 
That  is  the  body-induced  and  free-surface  induced  com¬ 
ponents  of  pressures  cancel  at  points  exterior  to  the 
hull  surface.  In  practice,  they  may  not  cancel  unless 
the  pressures  are  evaluated  numerically  in  certain  ways. 
It  has  been  found  that  using  the  free  surface  induced 
acceleration,  rather  than  velocity,  in  the  hull  normal 
boundary  condition  is  very  effective  in  removing  any 
influence  on  the  pressures  from  the  standing  waves 
generated  internal  to  the  hull. 

The  spectral  representation  of  the  free  surface  can  use 
either  a  rectangular  or  cylindrical  representation  in 
spectral  space.  That  is  n(x,  y),  the  elevation,  can  be 
represented  either  as 


NKX  NKY  ikx  • x  r 

n(x,y)  =  Y.  Akx  y  Aky  e  m  •  A 

n  m  L  nm 

n= l  m— i 


ikym-y 


+  A*  e 
nm 


-iky  -y 


m 


or 

n(x,y) 


NK  NTH 

y  k  -Ak  y  A9  A 
n  n  nm 

n=l  m— i 


ik  (x  cos0  +  y  sin0  ) 
n  m  '  m 


(01  2’  8NTH  +  2) 


Both  methods  have  been  programmed.  They  are  very  similar. 
The  advantage  of  the  rectangular  system  is  that  the 
length  scales  of  the  free  surface  are  independent  in  the 
x  and  y  directions.  This  is  particularly  useful  for 
typical  ship  hulls  where  the  hull  length  is  much  greater 
than  the  transverse  dimensions  of  the  ship  and  the  flow 
is  expected  to  change  much  more  slowly  in  the  x  direction 
than  in  the  y  direction. 
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Application  to  Linearized  Frequency-Domain  Seakeeping 

Characteristics  with  Forward  Speed 

Consider  a  hull  moving  with  forward  speed  U  and  oscillat¬ 
ing  about  its  mean  position  defined  by  the  surface 

F(x  -  Ut,  y,  z)  =  0 

The  flow  velocity  perturbation  induced  by  the  forward  motion  alone 
is  assumed  to  be  small,  i.e.  0(e).  Harmonic  perturbations  about 
the  mean  flow  field  are  also  assumed  to  be  of  order  e.  Interac¬ 
tions  between  the  steady  and  harmonic  portions  of  the  flow  field 
are  assumed  to  be  of  higher  order  in  e . 

Motions  are  described  in  a  coordinate  system  (x*,  y* ,  z*) 
moving  with  the  mean  hull,  i.e.  with  forward  speed  U.  The  coordi¬ 
nates  of  this  moving  system  relative  to  that  of  the  space-fixed 
(x,  y,  z)  system  are  given  by 

x*  =  x  -  Ut 

y*  =  y 

z*  =  z  . 


The  velocities  of  the  craft  relative  to  a  moving  (x*,  y* ,  z*)  ship- 
fixed  system  are  written  for  the  six  degrees  of  freedom  as 

v^  =  surge  velocity  (+  forward) 
v2  =  sway  velocity  (+  to  starboard) 
v3  =  heave  velocity  (  +  downward) 

v4  =  roll  velocity  (+•  for  right-hand  rotations  about  x  axis) 

Vg  =  pitch  velocity  (+  for  right-hand  rotations  about  y  axis) 

Vg  =  yaw  velocity  (  +  for  right-hand  rotations  about  z  axis) 

All  velocities  are  assumed  to  be  of  order  e  and  harmonic  in  time 
with  frequency  w.  [Since  these  velocities  are  specified  in  a 


ship-fixed  system,  they  are  related  to  both  velocities  and  angular 
displacements  in  a  space-fixed  system  for  finite  forward  speed. 

For  example, 


Vo(t)  =  V,(t)  -  0  c- (  t  )  •  u 

oHIP  SPACE 

where  U  is  the  forward  speed  and  0g  is  the  pitch  angle.] 

The  normal  velocity  at  a  point  x*  on  the  surface  is 

given  by  harmonic  velocities  in  a  ship-fixed  system  with  magnitude 

v  and  phase  0  is 
m  ^  m 


^  6 

v*  *  n  =  T,  n  (x*)  v  cosOt  +  0  ) 
mm  m 

m=l 

where  n^(x*)f  n^(x*),  and  n0(x*)  are  generalized  normals  defined  as 

n4(x*)  =  n2(x*)  •  x*  -  n3(x*)  •  x* 

•% 

n5(x*)  =  n3(x*)  •  x*  -  n1(x*)  •  x* 
n0(x*)  =  n1(x*)  •  x*  -  n2(x*)  •  x* 


similarly  the  normal  acceleration  produced  by  the  motion  is 

a*  •  n  =  -  >.  n  (x*)  v  u  sin(oit  +  0  ) 
mm  m 

m=l 

The  problem  is  to  compute  the  forces  and  moments  induced  by  these 
harmonic  motions  for  arbitrary  w  from  a  set  of  time-domain  solu¬ 
tions  (one  for  each  of  the  six  degrees  of  freedom)  for  impulsive 
motions.  Since  the  problem  is  linearized,  a  single  degree  of 
freedom  can  be  considered  without  loss  of  generality.  Consider, 
for  example,  the  steady  harmonic  motion  described  by 


n  =  v. 


nj(x*> 


coswt 


n  =  -<*>  v .  n.(x*)  simut 
J  J 


Assume  that  the  time-domain  solution  based  on  the  computational 
method  described  above  is  available  for  the  following  problem  in 
the  limiting  case  as  y  -*  0.  The  flow  is  at  rest  for  t  <  0.  Start- 
ing  at  time  t  =  0  the  flow  is  disturbed  by  motion  in  the  j  degree 

of  freedom  described  by  the  body  boundary  condition, 


n  =  n  .  ( x* ) 
J 


0  <  t 


=  0  t  <  0 


i.e.,  v*  •  n  =  n.(x*)  H(t) 

J 

where  H(t)  is  the  unit  step  function  and  where  x*  is  any  point 
on  the  hull  surface.  Define  the  resulting  potential  at  a  point 
x*  on  the  hull  as 

<(>.(x*,  t)  t  >  0 

«J 

and  the  forces  and  moments  relative  to  the  (x*,  y*,  z*)  system  in 
the  six  degrees  of  freedom  as 

FkJ(t)  k  -  1.  2,  3  ...  6  t  >  0 

The  tilda  symbol  ~  indicates  a  time-domain  solution  obtained  from 
the  numerical  simulation.  The  j  subscript  refers  to  the  degree  of 
freedom  of  the  motion.  The  k  subscript  refers  to  the  degree  of 
freedom  of  the  force  or  moment. 

The  forward  speed  U  should  be  included  in  the  time-domain 
computations.  Although  forward  speed  does  not  appear  directly  in 
the  body  boundary  condition,  it  enters  the  computation  in  two  ways. 
First  the  mean  hull  surface  moves 


F(x  -  Ut ,  y ,  z)  =  0 


with  speed  U  so  that  over  the  time  of  the  computation  the  space- 
fixed  (x,  y,  z)  and  body  fixed  (x*,  y*,  z*)  coordinate  are  dis¬ 
placed  with  respect  to  each  other.  Secondly,  the  body-induced 
pressure  contains  a  term  proportional  to  forward  speed- 


d  <f> 


PBD  =  _0 


BD 


34) 


dt 


+  pU 


BD 


3x 
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It  is  assumed  that  the  general  form  for  the  force  and 

moment  vector,  Fk(t)  generated  by  a  velocity  of  mode  j  and  time 

history  v.(t)  in  the  ship  fixed  reference  system, 

J 

v*  •  n  =  Vj(t)  nj (x)  , 


Fkj(t>  =  Fkj(IMP)  *<*>  +  /  Kkj^-T>  VT)  dT 

—  CO 

where  the  first  term  represents  virtual  moss  effects  and  the 
second  term  contains  the  memory  effects  resulting  from  the  free 
surface.  If,  for  example,  a  step  function  is  imposed  in  mode  j, 


v,  =  H(t) 

J 


then  the  solution  is 


t 

•/ 


Fkj(t)  =  /  K(r)dr  t  >  0 


and  Fkj(IMP)  is  defined  by 


\KIMP)  =  -P  ffiS B  V*'  °+)  nk(I)  ' 


From  integration  by  parts  for 


v*  •  n  =  v.  n.(x)  coswt 
3  3 


the  solution  is 


Fkj(t)  =  W  [Fkj  ( IMP)  Sinajt  +  J  Fkj  (T) 

0+ 


=  wsinwt  [Pkj(IMP)  +  J  Fk.(T)  cosiox  dx 

0+ 


-  wcoswt  J  Fkj(x)  sinwx  dT 


The  first  term  in  the  expression  for  Fkj  is  in  phase  with  the 
acceleration  and  may  be  thought  of  as  a  generalized  added  mass 
term.  The  second  term,  in  phase  with  velocity,  represents  a 
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damping.  The  hydrostatic  effect,  not  included  in  these  dynamic 
forces,  would,  of  course,  give  a  term  proportional  to  displacement. 


Wave  excited  forces  for  harmonic  waves  can  be  computed 
using  similar  techniques.  Consider  an  ambient  wave  field  defined 
by 

<t>  *  Re  I A  eilc^x  cos9  +  y  sin6)  -  kz  -  icot} 

He  |  a  eilc(x*  c°s9  +  y*  sinG  -  kz*  -  icoet  | 

where  A  is  the  complex  amplitude  and  coe  is  the  encounter  frequency 
defined  as 


w  =  oj  -  Uk  cosG  . 
e 

The  resulting  ambient  velocity  and  acceleration  field  can  be 
written  in  the  form 


and 


v(x*,  t)  =  v.(x*)  sin  to  t  +  v_(x*)  cos  w  t 

1  S  6 


a( x* ,  t)  =  ue  v^x*)  cos  coet  _  we  v2(x*)  sin  met 
=  a-^t)  vx(x*)  +  a2(t)  v2(x*) 


where 


and 


vx(x)  =  lm  |  A  ?  •  eik(x*  0030  +  y*  slne)  -  kz*  ! 


v2(x)  =  Re  I  A  V  •  eik(x*  cose  +  y*  slne)  '  kz*  I 


Now  assume  the  time-domain  method  of  solution  is  applied  to  the 
following  pair  of  time  step  problems.  In  each  case  the  hull  and 
the  surrounding  fluid  is  at  rest  for  t  <  0.  After  t  =  0;  the  body 
boundary  conditions  are 


v* 


•  n  =  vx(x*)  •  n(x*)  a^(t)  t  ^  0 


=  vx(x*)  n(x* )  H(t),  all  t 


A 
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for  the  first  case  and 


v*  •  n  =  v2(x*)  •  n(x*)  a2(t)  t  >_  0 
=  v2(x*)  •  n(x)  H(t),  all  t. 
i.e.,  ax(t)  =  a2(t)  =  H(t) 

for  the  second  case.  The  pair  of  time-domain  forces  and  moments 

for  each  of  the  six  degrees  of  freedom  obtained  by  applying  the 

time-domain  computational  solution  to  these  two  step  problems  for 

arbitrary  a, (t)  and  a0(t) 

1  *  t 

Fkl(t)  =  Fkl( IMP)  al(t)  +  f  Fkl(r)  a1(t“T)dT 

*X) 

t 

Fk2(t)  =  Fk2( IMP)  a2(t)  +  f  Fk2(T)  a2(t-T)dT 

'/0 

Where  the  pair  of  time-domain  forces  and  moments  for  each  of  the 
six  degrees  of  freedom  obtained  by  applying  the  time-domain 
computational  method  to  these  two  problems  are  written  as 

Fkl ( t )  and  ^(t) , 

respectively  where  the  tilda  ~  symbol  indicates  a  time-domain 
solution,  the  subscript  k  refers  to  the  degree  of  freedom  of  the 
force  or  moment,  and  the  remaining  subscript  indicates  which  of 
the  pair  of  body  boundary  conditions  were  employed. 

The  forces  and  moments  generated  by  the  harmonic  velocity 

field 


v(x*,  t)  = 


vx(x*) 


sin 


w  t 
e 


v2(x*) 


cos 


U)  t 

e 
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may  be  expressed  in  terms  of  the  time-domain  solutions  as 


Fk  =  «e/Jkl(T)  cos  oj0(t-T)  dr  -  .e/Fk2(T)  sin  we(t-T)  dT 
0  0 


F,  -  to  cosco^t 
k  e  e 


Put 


IMP) 


5  [^kl(T)  cos  weT  +  ^k2(r)  sin  a)eT]dT] 


0+ 

oo 


sinco  t  F 


e  L  k2( IMP) 


h 

0+ 


+  /^k2  (t)  COS  0JaT  -  Ki(T)  sin  oj0t 


kl 


]dT] 


where 


and 


?kl(IMP)  -Ovffte  tc)  nk(^) 

SB 

*k2(IMP)  *  *2<**'  ‘o’  V**) 
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Heave  Motion  of  a  Floating  Sphere 

As  a  simple  example  of  the  method  described  above,  the 
computer  program  capable  of  handling  any  three  dimensional  hull 
shape  was  applied  to  the  case  of  heave  motion  for  a  semi-submerged 
sphere,  i.e.,  a  hemispherical  hull.  The  computer  program  is  com¬ 
pletely  general  and  makes  no  assumptions  concerning  either  lateral 
or  longitudinal  symmetry  of  the  body  or  the  flow.  For  the  case 
of  zero  forward  speed,  however,  the  symmetry  of  the  hemispherical 
hull  requires  that  only  heave  force  result  from  heave  motion. 

The  magnitudes  of  other  forces  calculated  numerically  should  be 
one  measure  of  the  adequacy  of  the  surface  panels  chosen  to 
represent  the  body  and  the  program  itself. 

In  the  test  case  60  panels  were  used  to  represent  the 
hemisphere  with  spacings  specified  by 

Aip  -  2tt/12  and  A9  =  tr/5 

where  0  £  ^  <  2tt  is  the  circumferential  angle  and  0  £  0  £  it  is 
the  polar  angle.  The  hull-fixed  coordinates  are  specified  at 
each  of  the  four  corners  of  each  panel  with  coordinates  given  by 

x*  =  a  cos \p  cos0 
y*  =  a  sinip  cos9 

and  z*  -  a  sin9 

with  radius  a  set  to  unity  in  this  test  case. 

If  the  hull  is  not  flat  over  a  panel,  its  finite  size 
causes  errors  in  the  computation  of  the  panel  center  and  panel 
area.  For  simplicity  the  panel  center  is  computed  to  be  at  the 
mean  of  the  four  corner  points.  For  a  convex  hull  such  as  the 
hemisphere  these  points  fall  slightly  inside  the  physical  hull 
rather  than  directly  on  it.  More  important  is  the  error  in  the 
computed  hull  area.  Again  for  simplicity,  the  pair  of  diagonals, 

*31  *  x3  -  xx 
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and 


x42  =  x4  ‘  x2 

are  used  to  estimate  the  panel  area, 


and  unit  normal, 

n  =  x42  x  x31/|x42  x  x31| 

where  x^ ,  Xg,  x3  and  are  the  coordinates  of  the  four  corner 
points  of  the  panel  in  the  body-fixed  system.  Although  it  is 
immaterial  whether  any  single  panel  normal  points  inward  or  out¬ 
ward,  for  simplicity  all  normals  are  assumed  to  be  outward  in  the 
description  of  the  method  and  in  the  sample  problem.  In  the  case 
of  convex  bodies  the  areas  are  underpredicted  by  an  amount  which 
depends  on  the  panel  size.  Likewise,  the  areas  are  also  under¬ 
estimated  for  concave  portions  of  the  physical  hull. 

For  the  case  of  a  hemispherical  hull  represented  by  sixty 
plane  quadralateral  panels  the  impulsive  forces  and  moments  result¬ 
ing  from  an  instantaneous  change  in  the  heave  (z)  velocity 
(+  downward)  from  zero  at  t  =  0-  to  unity  at  t  =  0+,  was  computed 
to  be  (for  unit  density,  p  =  1): 

=  -0.00003 
=  .00656 

=  -1.00150 
=  . 00004 

=  -  .00011 
=  -  .00012  . 

From  theory 

Fk3(IMP)  =  0  k  *  3 

and 

F33(IMP)  *  -ff/3  *  -1' 04720  • 


F13(IMP)  =  surge  force 
F23( IMP)  =  sway  force 
F33( IMP)  =  heave  force 

F45( IMP)  =  ro11  moment 
F53( IMP)  =  Pitch  moment 

and  F63( IMP)  =  yaw  moment 
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The  computed  terms  are  thus  <  0(10  )  for  all  forces  and  moments 

other  than  heave.  The  computed  heave  force  is  about  3  percent 
too  low;  a  number  which  could  be  reduced  with  a  finer  hull  reso¬ 
lution  in  the  panels. 

The  time-domain  numerical  method  was  applied  to  the  semi- 
submerged  sphere.  Initially  (t  <  0)  the  flow  and  the  body  was 
specified  to  be  at  rest.  At  time  t  =  0,  a  step  function  was 
applied  to  the  heave  velocity.  Forward  speed  was  set  at  zero. 

Thus  the  hull  velocity  boundary  condition  is 


(vbd(x*)  +  vFS(x*))  •  n(x*)  =  nzU*)  H(t) 

where  x*  is  any  point  on  the  hull  surface  and  H(t)  is  the  unit 
step  function.  For  the  initial  impulse  time  step  of  infinitesmal 
duration  starting  at  t  =  0,  a  velocity  boundary  condition  was 
applied  at  each  panel  center  as 


j?!  Eij  °j(t  =  °+)  *  ^  •  “i  *  nzi 


(since  v^g  is  zero  at  t  =  0+).  The  source  strengths  at  t  =  0+, 
cK(t  =  0+)  were  then  used  to  find  the  potentials  on  the  panels 
and,  therefore,  the  forces  and  moments  acting  on  the  body  in  all 
six  degrees  of  freedom  as  a  result  of  the  initial  impulsive  step 
from  zero  to  unity  in  heave  velocity  (the  third  mode).  That  is 
equivalent  to 


Fkj(IMP)  k  “  l’  2’  ’ ' *  6  • 

In  practice  these  impulsive  forces  are  equivalent  to  the  formula 
given  earlier, 


k  j ( IMP ) 


=  -P  JJdS B  $j(x*,  0+)  nj 

ss 


(x*) 


For  later  time  steps  the  time-domain  force 


Fkj(t> 


I 


(T)dT,  t  >  0+ 
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was  computed  using  the  method  presented  in  this  report  with  an 
acceleration  body  boundary  condition  of  the  form 


nb 

jS  Eij  0j(t>  =  “*(t)  -  V” 

(with  the  forward  speed  U  set  to  be  zero) 


i 


TT  V 

3x  Vxs 


t  >  0+ 


since  the  prescribed  normal  acceleration  a|  is  identically  zero 
for  t  >  0+.  The  resulting  values  for  o^( t),  i.e.,  the  time  deriva 
tive  of  the  source  strength  density  on  each  panel,  gives  the  time 
rate  of  change  of  the  body  potential  and,  therefore,  the  body- 
induced  pressure  Pg^.  at  each  panel.  By  adding  the  free-surface 
induced  pressures,  the  total  pressure  acting  on  each  panel  was 
determined.  This  produced  the  net  computed  force 


(t) 


acting  on  the  body  for  t  >  0+. 


For  the  numerical  test  case,  the  forces  and  moments  of 
the  class  F^Ct),  k  f  3  were  all  insignificant.  Figure  1  shows 
F33(t)  computed  under  the  following  numerical  conditions  - 


L  =  maximum  length  scale  in  x  =  2.5 

Z  =  minimum  length  scale  in  x  -  0.25 

Ly  =  maximum  length  scale  in  y  *  2.5 

Ay  *  minimum  length  scale  in  y  =  0.25 

T  =  maximum  time  scale  =  8.0 

g  =  acceleration  of  gravity  =  1.0 

p  =»  fluid  density  *1.0 

a  *  radius  of  sphere  *  1.0 

At  *  time  step  size  *  0.  10 
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rORCE 


FORCE  HISTORY 


Figure  1.  Time  history  of  force  history  induced 
by  a  unit  step  in  heave  velocity  start 
ing  at  t“0  for  a  semisubmerged  sphere. 
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The  first  five  of  these  conditions  produce  17  wave  lengths  for 
the  free  surface  representation  in  both  the  x  and  y  directions  or 
17  x  17  modes.  For  the  harmonic  case 


Vg(t)  =  coscut  , 


the  force  acting  in  the  z  direction  (+  downward)  is  given  by 
Fg( t )  =  -ui  sin(uit)  f  F, 


00 

33(  IMP)  +J  F33(T)  cos(u)T)dT 
0+ 


OC 

/ 


+ui  cos(uit)  J  Fgg(-r)  sin(cuT)dT  . 

0+ 

This  formula  can  be  evaluated  numerically  and  compared  with  the 
classical  results  of  Havelock  [3] . 

Havelock  expresses  the  z  force  induced  by  the  heave 
oscillations 


v3  *  cosuit 


a-s 

23  2  3 

Z  =  g-  irpa  k  sinwt  -  -g  irpa  2h  ui  cosuit 

where  a  is  the  radius  and  p  is  the  density  (both  are  unity).  The 
coefficients  k  (not  a  wave  number)  and  2h  determine  the  added  mass 
and  hydrodynamic  damping  of  the  sphere  in  heave.  These  two  coef¬ 
ficients  are  commonly  expressed  in  terms  of  the  nondimensional 
2 

variable  ui  a/g. 

These  coefficients  defined  by  Havelock  are  related  to 
the  values  derived  from  the  time-domain  solution  for  a  step  func¬ 
tion  in  the  heave  mode  as 

00 

k(ui  a/g)  2rf  F33(  IMP)  +  f  F33^  cos(WT)dT 

0+ 
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and 

oo 

2h(oj2a/g)  =  -  ^  y*  Fgg( t )  sin(taT)dT 

0+ 

2  2 

Figure  2  shows  coefficients  k(<o  a/g)  and  2h((o  a/g)  calculated  in 

this  manner  using  ~  _1-0015  and  F33(t)  as  shown  an  Figure  1. 

The  resulting  curves,  close  to  those  predicted  by  Havelock  using 

analytic  methods,  given  in  reference  [3] .  Tabular  results  for 

this  problem  were  also  given  by  Kim  [4] .  Table  1  below  shows  a 

comparison  between  Kim's  results  and  the  computed  results  for  var- 

2 

ious  values  of  to  a/g. 
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Table  1.  Comparison  of  analytic  and  numerical 
results  for  added  mass  and  damping 
coefficients  of  a  semisubmerged  sphere 
of  radius  a  oscillating  with  frequency  w . 


Added  Mass  Coefficients 


2  / 
a/g 

k(ai2a/g)  [Ref.  [4]] 

2 

k(ou  a/g)  [Computed] 

0.50 

0.60 

0.58 

1.00 

0.43 

0.43 

1.50 

0.39 

0.39 

2.00 

0.39 

0.39 

2.50 

0.40 

0.39 

3.00 

0.41 

0.38 

Damping  Coefficients 


2 

a)  a/g 

2h(w2a/g)  [Ref.  [4]] 

2 

2h(w  a/g)  [Computed] 

0.50 

0.34 

0.33 

1.00 

0.25 

0.25 

1.50 

0.17 

0.17 

2.00 

0.11 

0.12 

2.50 

0.07 

0.10 

3.00 

0.06 

0.06 
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MASS  AND  DAMPING  COEFFICIENTS 
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